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Hodographic  analysis  of  gravity  waves: 

Relationships  among  Stokes  parameters,  rotary  spectra 
and  cross-spectral  methods 

Stephen  D.  Eckermann1 

Computational  Physics,  Incorporated,  Fairfax,  Virginia 

Abstract.  The  elliptical  rotation  with  height  of  horizontal  velocities  produced  by  gravity 
waves  provides  considerable  information  about  the  wave  field.  Methods  of  statistically 
characterizing  velocity  ellipses  from  data  currently  fall  into  three  main  categories:  (1) 
hodographic  analyses,  (2)  cross-spectral  analyses,  and  (3)  rotary  spectral  analyses.  The 
three  methods  have  some  intuitive  similarities,  yet  precise  interrelationships  among  them 
are  presently  unclear.  The  three  techniques  are  interrelated  here  using  the  so-called 
“Stokes  parameters”  of  the  wave  field,  which  initially  provide  a  concise  description  of  the 
hodographic  analysis  method  (1).  On  Fourier  transforming  the  Stokes  parameters, 
standard  formulae  employed  in  rotary-spectral  and  cross-spectral  analysis  methods  can 
then  be  expressed  in  terms  of  the  resulting  “Stokes-parametcr  spectra.”  The  results 
highlight  some  drawbacks  in  the  use  of  cross-coherence  spectra  between  velocity 
components  to  verily  the  existence  of  a  coherent  wave  motion.  A  more  robust  measure  is 
suggested,  based  on  the  “degree  of  polarization”  of  classical  hodograph-based  Stokes- 
parameter  analysis,  which  can  be  generalized  to  provide  an  analogous  spectral  measure  for 
evaluation  in  rotary-spectral  and  cross-spectral  analyses. 

Introduction 

It  is  well  known  that  the  horizontal-velocity  motion  ellipses 
produced  by  gravity  waves  can  provide  a  lot  of  information  on 
the  wave  field.  For  a  single,  linear,  steady,  nondissipating, 
plane  wave  in  a  region  with  negligible  vertical  shear  in  the 
background  wind,  the  axial  ratio  of  this  motion  ellipse  can 
provide  the  intrinsic  frequency  of  the  wave,  a  parameter  which 
is  very  difficult  to  measure  any  other  way  [e.g.,  Cot  and  Barai, 

1986] .  The  horizontal  propagation  axis  of  the  wave  is  given  by 
the  orientation  of  the  ellipse  major  axis  [Vincent  and  Fritts, 

1987] ,  while  the  sense  of  ellipse  rotation  with  height  provides 
the  direction  of  the  wave’s  vertical  group  velocity  [e.g.,  Hirota 
and  Niki,  1985;  Andrews  el  at.,  1987],  In  strongly  sheared  en¬ 
vironments,  ellipticity  may  also  arise  due  to  the  interaction  of 
the  wave  motion  with  the  background  velocity  profile  [Hines, 

1989;  Cho,  1995].  Consequently,  a  number  of  techniques  have 
been  developed  which  can  quantify  the  elliptical  properties  of 
vertical  profiles  of  horizontal  velocity  perturbations  produced 
by  gravity  waves. 

The  first  and  most  widely  used  method  relics  on  forming 
hodographs  from  the  velocity  profiles.  This  is  a  straightforward 
and  accurate  technique  when  a  large-amplitude  monochro¬ 
matic  wave  is  present  in  a  region  of  minimal  background  wind 
shear  [e.g.,  Cot  and  Burnt,  1 986;  Muraoka  et  at.,  1994],  How¬ 
ever,  actual  velocity  hodographs  are  usually  irregular  and  vari¬ 
able,  and  so  a  number  of  hodographic  algorithms  have  been 
developed  which  sort  and  classify  a  large  number  of  profiles  to 
yield  a  stalistical  picture  of  Ihe  wave  field  [e.g..  Hirota  and  Niki, 
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1985;  Hamilton,  1991;  Hall  et  al.,  1995].  One  particularly  gen¬ 
eral  form  of  the  hodographic  method  is  a  description  in  terms 
of  Stokes  parameters  [Vincent  and  Fritts,  1987],  an  approach 
we  shall  pursue  here. 

However,  observations  and  models  have  suggested  that  a 
fairly  broad  spectrum  of  many  waves  is  present  on  average 
[e.g.,  Tsnda  et  al.,  1991;  Hines,  1991],  whereas  hodographic 
techniques  require  a  single  coherent  wave  to  dominate  the 
fluctuations  within  the  profile.  Eckermann  and  Hocking  [1989] 
demonstrated  that  hodographic  analyses  of  profiles  which  con¬ 
tained  a  large  number  of  gravity  waves  sometimes  gave  results 
which  reflected  the  statistical  properties  of  random  incoherent 
fluctuations,  rather  than  the  polarization  characteristics  of  any 
waves  in  the  data.  Consequently,  spectral  techniques  have  been 
pursued  which  can  search  for  coherent  waves  within  various 
wavelength  bands  and,  if  any  are  identified,  can  isolate  each 
wave  and  characterize  it  using  hodographic  methods. 

First  efforts  along  these  lines  concentrated  on  rotary  spectra, 
which  are  power  spectra  of  the  complex  velocity  vector  u'(z ) 
+  iv'(z),  where  u'(z)  and  v'(z)  are  vertical  profiles  of  the 
fluctuating  zonal  and  meridional  wind  velocities,  respectively. 
When  plotted  in  the  complex  plane,  the  rotation  of  this  “rotary 
vector”  with  height  provides  identical  information  to  velocity 
hodographs.  Thus  generalized  rotary-vector  analysis  methods 
have  been  developed  [e.g.,  Gonellu,  1972;  Hayashi,  1979]  and 
are  used  commonly  in  studies  of  oceanic  gravity  waves  [e.g., 
Artale  and  Gasparino,  1990;  Plueddemann,  1992].  When  ap¬ 
plied  to  gravity  waves  in  the  middle  atmosphere,  strong 
anisotropies  between  the  clockwise-rotating  and  anticlockwise- 
rotating  rotary  spectral  components  were  revealed  at  all  verti¬ 
cal  wavenumbers,  which  were  consistent  with  a  preponderance 
of  upward  propagating  gravity  wave  packets  [e.g.,  Thompson, 
1978;  Cadet  and  Teitelhaum,  1979;  Vincent,  1984;  Hass  and 
Meyer,  1987;  Eckermann  and  Vincent,  1989;  Thomas  et  al., 
1992], 
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ECKERMANN:  RELATIONSHIPS  AMONG  GRAVITY  WAVE  ANALYSIS  METHODS 


Eckermann  and  Vincent  [1989]  presented  a  spectral  version 
of  the  Stokes-parameters  method  of  characterizing  gravity 
wave  hodographs  [Vincent  and  Fritts ,  1987]  and  applied  it  to 
rocketsonde  profiles  of  stratospheric  velocities.  The  technique 
provided  a  “degree  of  polarization”  measure  for  identifying 
coherent  waves  within  any  given  wavenumber  band  as  well  as 
other  hodograph-related  measures  of  the  polarization  charac¬ 
teristics  of  coherent  waves  within  various  bands  [Eckermann 
and  Vincent,  1989;  Eckermann  el  al.,  1995], 

Recently,  Cho  [1995]  has  developed  a  cross-spectral  tech¬ 
nique,  which,  he  argued,  combined  the  advantages  of  the  mul¬ 
tiwave  rotary-spectrum  approach  with  the  interpretive  advan¬ 
tages  of  a  hodographic  analysis.  On  the  latter  property  he 
noted  that  the  technique  was,  in  principle,  a  spectral  analogue 
of  the  Stokes-parameters  analysis.  He  employed  the  method  to 
investigate  the  elliptical  properties  of  the  wave  field  measured 
by  radar  over  Arecibo. 

Clearly,  all  three  spectral  methods  give  similar  results  and 
must  be  related  in  some  way.  The  purpose  of  this  short  paper 
is  to  explicate  these  relationships. 


Hodographic  Decomposition:  Stokes  Parameters 

In  hodographic  analyses  of  gravity  wave  horizontal-velocity 
perturbations,  velocities  are  plotted  vectorally,  and  values  at 
successive  heights  are  interconnected.  The  resulting  velocity 
hodographs  frequently  reveal  elliptical  wave-induced  velocity 
“spirals”  as  a  function  of  height  [Hirota  and  Niki,  1 985 ;  Cot  and 
Barat,  1986;  Yamanaka  et  al.,  1996].  Various  analysis  tech¬ 
niques  have  been  developed  to  identify  and  characterize  large 
numbers  of  these  hodographs  in  some  statistical  sense  [ Hirota 
and  Niki,  1985;  Kitamura  and  Hirota,  1989;  Hamilton,  1991; 
Hall  et  al.,  1995], 

One  complete  and  concise  formulation  of  this  method  is 
provided  by  the  so-called  Stokes  parameters.  This  is  a  standard 
analysis  technique  for  electromagnetic  waves  and  is  well  de¬ 
scribed  in  standard  texts  from  optics  and  radiophysics  [e.g., 
Kraus,  1966;  Bom  and  Wolf. ,  1980].  An  attractive  feature  of 
these  parameters  is  that  any  partially  polarized  wave  motion 
can  be  described  completely  and  uniquely  in  terms  of  its  Stokes 
parameters. 

Vincent  and  Fritts  [1987]  were  the  first  to  advocate  their 
usefulness  in  gravity  wave  studies.  In  this  context  the  technique 
assumes  that  any  given  vertical  profile  of  horizontal-velocity 
perturbations,  («'(z),  t/(z)),  contains  a  partially  polarized 
wave  field:  that  is,  a  single  coherent  wave  is  present,  of  peak 
amplitude  («n,  u0),  within  an  unpolarized,  isotropic,  noiselike 
background  velocity  field  of  variance  The  four 

Stokes  parameters  can  then  be  defined  as 

f  =  \(u})  +  »o)  +  («Lse  +  Vnoise)  =  u'2  +  v'Z,  (1) 

=  =  (2) 

P  =  k,)U,)  cos  8  =  2 u'v' ,  (3) 

Q  =  uaV( ,  sin  5,  (4) 

where  overbars  denote  height  and/or  time  averaging  of  the 
perturbation  velocities;  and  S  is  a  phase  defining  the  ellipticity 
of  the  wave:  8  =  0°  or  180°  implies  linear  polarization  (Q  = 
0),  8  —  90°  or  270°  implies  circular  polarization  (P  —  0),  and 
anything  in  between  implies  elliptical  polarization.  Thus  P  is 
the  “in-phase”  covariance  associated  with  linear  polarization, 


and  Q  is  the  “in  quadrature”  covariance  associated  with  circu¬ 
lar  wave  polarization,  while  I  clearly  quantifies  the  total  vari¬ 
ance  and  D  defines  its  axial  anisotropy. 

Note  that  the  isotropic  noiselike  variance  in  the  profile  does 
not  contribute  to  D,  P,  or  Q  but  does  contribute  to  the  total 
variance  I.  To  discriminate  between  the  wave  and  the  noiselike 
contributions  to  /,  a  degree  of  polarization  d  can  be  defined 
which  quantifies  the  fractional  contribution  of  the  coherent 
wave  motion  to  the  total  velocity  variance.  Since  I2  —  D2  + 
P2  +  Q2  for  a  single  wave  in  the  absence  of  noise,  then  it 
follows  that 


d  = 


(D2  +  P2+  Q2)'12 
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(5) 


where  Osrfsl  and  d  is  rotationally  invariant,  in  that  its 
value  does  not  depend  on  the  coordinate  axes  used  to  evaluate 
(l)-(4).  A  number  of  other  useful  parameters  can  also  be 
calculated  from  the  Stokes  parameters  [e.g.,  Kraus,  1966;  Bom 
and  Wolf,  1980;  Eckermann  and  I'incent,  1989], 

Most  Stokes-parameter  definitions  are  a  factor  of  2  larger 
than  those  in  (l)-(4)  [Vincent  and  Fritts,  1987],  The  inclusion  of 
some  constant  scaling  factor  has  no  influence  on  the  technique, 
and  we  adopt  the  present  definitions  because  they  make  later 
correspondences  with  spectral  quantities  particularly  simple. 


Stokes-Parameter  Spectra 

The  expressions  (l)-(4),  as  given  by  Vincent  and  Fritts 
[1987],  involve  computing  variances  and  covariances  front  the 
velocity  profile.  This  has  some  disadvantages:  for  example, 
computing  Q  from  definition  (4)  is  not  straightforward  and  was 
not  attempted  by  Vincent  and  Fritts  [1987]  ( Eckermann  et  al. 
[1995J  noted  that  it  can  be  achieved  by  Hilbert  transforming 
one  of  the  profiles).  This  profile-based  form  of  the  technique 
also  assumes  that  one  dominant  wave  exists  in  the  data, 
whereas  in  atmospheric  profiles,  it  is  more  usual  to  encounter 
a  number  of  superposed  waves  [e.g.,  Eckermann  and  Hocking, 
1989;  Riggin  et  al.,  1995]. 

Analysis  of  profiles  which  contain  more  than  one  wave  is 
clearly  more  suited  to  the  Fourier  domain,  and  Vincent  and 
Fritts  [1987]  noted  that  Q  as  well  as  the  other  parameters  might 
be  calculable  spectrally.  This  possibility  was  investigated  and 
developed  by  Eckermann  and  Vincent  [1989].  Given  vertical 
profiles  u'(z)  and  v’(z),  as  before,  then  Fourier  transforming 
them  over  their  full  height  ranges  yields  the  following  Fourier 
representations  of  the  profiles: 

U  (m)  =  UR(m)  +  tUi(m),  (6) 

V(m)  =  VR(m)  +  iVi(m),  (7) 

where  m  is  the  vertical  wavenumber,  and  the  bold  symbols 
denote  a  complex  term.  From  here,  Eckermann  and  Vincent 
[1989]  derived  the  following  power  spectral  densities  for  the 
four  Stokes  parameters,  based  on  definitions  (l)-(4): 

l(m)  =  A  (UR(m)  +  Uj(m)  +  V\{m)  +  Vj{m)),  (8) 

Dim)  =  A(U\{m)  +  Vj(m)  -  V2K(m)  -  F/(m)),  (9) 

P(m)  =  2  A(UR{m)VR(m)  +  U,{m)V,(m)),  (10) 

Q(m)  =  2 A{UR{,n)V,(m)  -  U,(m)VR(m)).  (11) 
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The  overbars  denote  averages  over  some  number  of  indepen¬ 
dent  spectral  realizations,  to  remove  the  effects  of  incoherent 
motions  (for  some  types  of  data,  this  is  either  inappropriate  or 
impossible),  and  A  is  a  constant  which  scales  the  squared 
Fourier  terms  to  power  spectral  densities.  Consequently,  any 
Stokes  parameter  X  (one  of  /,  D ,  P  or  Q)  can  now  be  eval¬ 
uated  over  the  full  range  of  heights  within  any  wavenumber 
band,  according  to  Parseval’s  relation 


=  ‘ X(m)dm ,  (12) 

mi 

where  Xmi  mi  is  one  of  the  Stokes  parameters  (1  )~(4)  evalu¬ 
ated  over  a  wavenumber  band  between  m,  and  m2  from  its 
corresponding  power  spectral  density  X(m).  This  Fourier  de¬ 
composition  of  the  Stokes  parameters  clearly  allows  for  more 
than  one  coherent  wave  to  be  isolated  from  any  given  sequence 
of  velocity  profiles,  so  long  as  the  various  coherent  waves  in  the 
profile  exist  within  well-defined  and  well-separated  vertical- 
wavenumber  bands. 

Note  that  (8)  is  just  the  total  horizontal-velocity  autospec¬ 
trum,  (9)  is  the  difference  between  the  component  autospectra, 
and  (10)  and  (I  l)  are  twice  the  real  and  imaginary  components 
of  the  total  velocity  cospectrum  (as  discussed  below).  These 
spectral  expressions  are  closely  related  to  the  various  compo¬ 
nents  of  the  classical  “coherency  matrix,”  which  can  also  be 
used  to  quantify  the  Stokes  parameters  [ Bom  and  Wolf,  1980]: 
indeed,  a  direct  analogy  was  provided  by  Caiman  [1978a] 
through  a  “spectrum  density  matrix,”  which  he  used  to  char¬ 
acterize  various  types  of  polarized  oceanic  wave  motions.  The 
key  point  here  is  that  when  coherent  waves  are  present  and 
well  separated  from  each  other  in  Fourier  space,  standard 
autospectra  and  cospectra  of  the  profiles  provide  fundamental 
information  about  the  polarization  characteristics  of  these 
waves.  The  usefulness  of  this  analogy  is  demonstrated  in  the 
next  two  sections,  where  some  other  spectral  formulae  which 
are  commonly  used  to  analyze  gravity  wave  fields  are  expressed 
in  terms  of  the  “Stokes-parameter  spectra”  (8)— (1 1). 


Stokes-Parameter  Spectra 
and  the  Cross-Spectral  Method 

Of  particular  interest  here  are  the  expressions  for  P(m )  and 
Q(m ),  which  are,  respectively,  twice  the  “in-phase”  and  “in- 
quadrature”  correlation  spectra  between  the  vertical  profiles 
w'(z)  and  u'(z).  Because  of  their  orthogonality  implied  by  (3) 
and  (4),  they  can  be  combined  into  a  complex  vector 

C(m)  =  P{m)  +  iQ(m)  -  T{m)  exp  [t€>(m)],  (13) 

where  C(m)  is  twice  the  total  cross  spectrum  between  u'(z) 
and  v'(z),  T2(m)  —  P2(m)  +  Q2(m),  and 

/  Q(m)\ 

<I>(w)  =  arctan  I  j.  (14) 

We  then  normalize  (13)  to  form  the  cross-coherence  spectrum 
[e.g.,  Fofonoff,  1969] 


pressed  entirely  in  terms  of  the  Stokes-parameter  spectra, 
yielding 


S(m)  = 


/  P\m)  +  Q\m )  \ 
\  -  D2(m)  ) 


1/2 

exp  [t<I>(m)]. 


(16) 


Thus  the  correspondence  between  Cho’s  cross-spectral 
method  and  the  hodograph-based  Stokes-parameters  method 
is  given  by  (14)  and  (16).  Several  points  are  now  worth  noting. 

First,  the  cross-spectral  method  of  Cho  [1995]  makes  most 
use  of  the  cross-spectral  phase  <t>(m).  Note  that  it  is  the 
spectral  analogue  of  the  S  term  in  (3)  and  (4),  as  Eckermann 
and  Vincent  [1989]  demonstrated  (see  their  equation  (7)).  On 
assuming  —180°  5  $(m)  180°,  then  Cho  [1995]  argued 

that  the  sign  of  <3>(m)  yielded  the  sense  of  elliptical  velocity- 
vector  rotation  with  height,  whereas  in  Stokes-parameter  anal¬ 
yses,  the  sign  of  Q(m)  is  used  [ Eckermann  and  Vincent,  1989; 
Eckermann  et  al.,  1995].  The  two  approaches  are  seen  to  be 
equivalent  on  inspecting  (14),  since  the  sign  of  4>(m)  is  deter¬ 
mined  solely  by  the  sign  of  Q(m). 

Cho  [1995]  also  advocated  using  the  magnitude  of  S(m ) 
(where  S(m)  has  been  suitably  averaged  over  a  number  of 
independent  spectral  realizations)  to  verify  the  presence  of  a 
coherent  wave  event  within  any  spectral  band.  This  requires 
both  height  and  time  resolution,  so  that  spectra  calculated 
from  height  (time)  profiles  can  be  averaged  with  respect  to 
time  (height).  It  is  more  usual  to  investigate  the  cross  coher¬ 
ence  between  two  time  series  of  a  given  quantity  which  arc 
recorded  at  displaced  locations,  to  infer  a  “coherence  length” 
of  the  recorded  phenomenon  [e.g.,  Madden  and  Julian,  1972; 
Briscoe,  1975;  Pinkel,  1975;  Reid  and  Vincent,  1987], 

|S(m)|,  on  the  other  hand,  quantities  the  coherence  be¬ 
tween  the  two  orthogonal  components  of  a  single-vector  time 
series.  This  leads  to  problems  which  are  analogous  to  those 
encountered  when  one  uses  the  correlation  coefficient  between 
w'(z)  and  v'(z )  (a  normalized  version  of  P  in  (3))  to  infer 
azimuthal  wave  field  anisotropy,  as  discussed  in  section  3.2  of 
Eckermann  and  Vincent  [1989],  As  an  illustration  of  the  prob¬ 
lem,  consider  a  profile  containing  a  single  linearly  polarized 
wave  aligned  either  zonally  or  meridionally,  which  implies  P  = 
(2  =  0.  According  to  (13)  this  wave  has  a  vanishing  cross 
spectrum  and  therefore  an  undefined  cross  coherence  (15) 
(adding  some  noise  eliminates  this  singularity,  giving  a  vanish¬ 
ing  cross  coherence).  Yet  clearly,  a  100%-coherent  wave  mo¬ 
tion  is  present.  This  weakness  in  the  cross-coherence  technique 
can  be  avoided  by  rotating  the  coordinate  axes  through  suc¬ 
cessive  angles  and  recalculating  the  cross-coherence  spectra 
each  time  [Eckermann  and  Vincent,  1989]. 

To  illustrate  the  problem,  Figure  1  compares  the  cross- 
spectral  and  Stokes-parameter  methods  of  inferring  coher¬ 
ence,  based  on  synthesized  profiles  containing  a  coherent  wave 
plus  noise.  In  Stokes-parameter  analysis  a  coherent  wave  is 
identified  through  the  degree  of  polarization  (5),  so  that, 
within  a  wavenumber  band  ml  s  m  £  m2, 


d,„im  = 


(Dl 


Q 


2  \  1/2 


(17) 


S(m)  =  2A\C(m)\  |V(m)|  '  (15) 

S(/n)  is  used  in  the  cross-spectral  analysis  method  of  Cho 
[1995],  With  the  aid  of  (8)— (11)  and  (13),  (15)  can  be  reex- 


Note  from  Figure  1  that  both  techniques  give  the  same 
results  when  the  semimajor  axis  of  the  wave  ellipse  is  aligned 
to  45°  to  the  coordinate  axes.  However,  when  the  same  wave  is 
oriented  either  north-south  or  east-west,  the  cross-coherence 
method  loses  sensitivity,  whereas  dm  ^  remains  insensitive  to 
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Measures  of  Wave-Variance  Content  versus  Azimuth 


Figure  1.  Comparison  of  the  cross  coherence 
=  S(m)  dm\ )  and  degree  of  polarization 
as  computed  from  synthesized  data.  The  synthesized  vertical 
profiles  contained  an  elliptic-ally  polarized  wave  (|//«|  =  0.1, 
where  (o  is  intrinsic  frequency  and /  is  inertial  frequency)  with 
a  total  horizontal-velocity  variance  of  25.25  m2  s  2  and  a  ver¬ 
tical  wavelength  of  4  km.  The  semimajor  axis  of  this  coherent 
wave  was  aligned  at  an  azimuth  angle  d>  and  then  sampled  over 
a  total  height  range  of  40  km,  and  at  a  vertical  resolution  is  1 
km  to  give  two  profiles  of  zonal  and  meridional  velocities. 
Uneorrelated  Gaussian  noise  with  a  variance  of  10  m2  s  2  was 
added  to  both  the  zonal  and  the  meridional  profiles.  One 
hundred  independent  pairs  of  zonal  and  meridional  velocity 
profiles  were  analyzed  spectrally  within  a  wavenumber  band 
between  m,  =  2w  (5  kni)  1  and  m2  =  2ir  (3  km)-1. 


the  wave’s  azimuthal  alignment.  These  properties  of  |S(m) 
and  dm  can  be  used  to  relate  the  two,  as  noted  in  Appen¬ 
dices  B  and  C  of  Hayashi  [1979],  The  relation  in  terms  of 
Stokes-parameter  spectra  is 


|S(«i)|2  = 


I2(m)d2(m)  —  D2(m) 
I2(m)  —  D2{m) 


(18) 


where 


-  D2(m)  +  P2(m)  +  Q2(m) 

rf(m)  = - TV) - ’  (t9) 

Note  that  (18)  reproduces  the  earlier  result:  that  a  linearly 
polarized  wave  aligned  either  north-south  and  east-west 
( P(m )  =  Q(m)  —  0)  has  an  undefined  cross  coherence  in  the 
absence  of  noise  ( d(m )  =  1)  and  a  vanishing  cross  coherence 
when  some  noise  variance  is  present  (0  <  d(m )  <  1). 

The  sensitivity  of  |S(/ti  )  |  to  ellipse  orientation  was  high¬ 
lighted  in  oceanographic  studies  by  Fofonoff  [1969],  who  pre¬ 
sented  some  alternative  measures  which  were  invariant  under 
coordinate  rotations.  This  stimulated  efforts  to  identify  other 
quantities  which  were  similarly  insensitive  to  coordinate  spe¬ 
cifics  and  led  to  the  development  of  the  rotary-vector  method 
[Gonella,  1972],  which  we  investigate  next. 


Stokes  Parameters  and  the  Rotary  Spectrum 

The  Stokes  parameters  are,  of  course,  not  the  only  way  in 
which  gravity  wave  velocity  hodographs  can  be  objectively 
characterized  [e.g.,  Hamilton ,  1991],  Another  quite  general 


approach  is  the  rotary-vector  method,  which  has  proved  a 
popular  tool  in  oceanic  studies  of  gravity  waves  [&.g.,Artale  and 
Gasparino,  1990;  D’Asaro  and  Morehead ,  1991;  Plueddemann, 
1992],  The  rotary-vector  method  utilizes  the  complex  velocity 
phasor  u'(z )  +  iv'(z)  and  provides  a  number  of  different 
parameters  which  describe  various  aspects  of  the  wave  polar¬ 
ization  [e.g.,  Gonella,  1972;  Caiman,  1 978b].  However,  in  at¬ 
mospheric  studies  to  date,  only  the  rotary  spectrum  has  been 
calculated  with  any  sort  of  regularity  [e.g.,  Thompson,  1978; 
Cadet  and  Teitelbaum,  1979;  Vincent,  1984;  Hass  and  Meyer, 
1987;  Thomas  et  at.,  1992],  and  so  we  shall  focus  mostly  on  this 
quantity.  Treatises  on  the  rotary-vector  method  and  its  corre¬ 
spondences  with  various  hodographic  and  cross-spectral  terms 
are  given  by  Gonella  [1972],  Mooers  [1973],  Caiman,  [1978a], 
and  Hayashi  [1979], 

Rotary  spectra  are  derived  by  first  Fourier  transforming 
u'(z )  +  lv’(z),  yielding  a  Fourier  transform  R(m).  Clock¬ 
wise-rotating  and  anticlockwise-rotating  circular  phasors  then 
partition  to  the  negative  and  positive  wavenumbers,  respec¬ 
tively,  of  R(m),  so  that  one  can  evaluate  clockwise-rotating 
(CW(m))  and  anticlockwise-rotating  (ACW(m))  power  spec- 


tral  densities  according  to 

A 

CfV(m)  =  y  R*(—m)R(—m), 

(20) 

-  A 

ACW (m)  =  y  R*(w)R(m), 

(21) 

where  asterisks  denote  complex  conjugation  [Leaman  and  San¬ 
ford,  1975;  Vincent,  1984],  The  factor  of  two  in  (20)  and  (21)  is 
needed  since  both  positive  and  negative  harmonics  of  the  Fou¬ 
rier  transform  are  used  here,  whereas  only  one  side  of  the 
Fourier  transform  is  used  in  (8)— (11)  because  (6)-(7)  are  Her- 
mitian.  The  total  autospcctrum  (equal  to  l(m))  is  given  by 
CW(m)  +  ACW(w)  [ Leaman  and  Sanford,  1975]. 

Standard  formulae  interrelate  the  amplitudes  of  the  zonal 
and  meridional  time  series  to  clockwise-rotating  and  anticlock¬ 
wise-rotating  circular  velocity  phasors,  and  Fourier  transform¬ 
ing  them  yields  relations  among  the  rotary  spectral  compo¬ 
nents  and  the  autospectra  and  cospectra  of  u'(z)  and  v'(z) 
[Gonella,  1972;  Mooers,  1973;  Caiman,  1978a].  Given  the  def¬ 
initions  (8)— (11)  of  the  Stokes-parameter  spectra  and  their 
earlier  stated  connections  with  autospectra  and  cospectra,  then 
these  results  can  be  used  to  relate  the  rotary  spectra  to  the 
Stokes-parameter  spectra.  The  results  are 

CW{m)  =  \U{m)-  Q(m)],  (22) 

ACW(m)  =  \l(m)  +  Q(m)).  (23) 

Note  that  l(m)  =  CW(m )  +  ACW(m),  as  required. 

The  above  expressions  have  some  practical  advantages.  For 
example,  a  conventional  rotary  spectrum  estimate  involves  fil¬ 
tering  the  raw  data  using  fast  Fourier  transform  (FFT)  meth¬ 
ods  (two  FFTs,  followed  by  two  inverse  FFTs)  and  then  re¬ 
transforming  u'(z)  +  lv'(z)  [e.g.,  Vincent,  1984],  a  total  of 
five  FFTs.  A  rotary  spectrum  calculated  directly  from  the 
Stokes-parameter  spectra  (8)  and  (11)  of  the  raw  data,  using 
(22)  and  (23),  requires  only  two  FFTs. 

As  in  the  cross-spectral  method,  it  is  also  common  in  the 
rotary-vector  method  to  evaluate  a  cross-coherence  spectrum 
S±(«i)  between  clockwise-rotating  and  anticlockwise-rotating 
velocity  components  [e.g.,  Fofonoff,  1969;  Gonella,  1972;  Ha- 
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yashi,  1979].  This  coherence  measure  can  also  be  expressed 
solely  in  terms  of  Stokes-parameter  spectra,  as  follows: 

/  D\m)  +  P2(m)\  1/2 

=  2V))  “p  (2t) 

where 


4,=(m)  =  arctan^j.  (25) 

Analogous  derivations  in  terms  of  autospectra  and  cospectra 
are  given  by  Gonella  [1972]  and  Mooers  [1973].  Similar  rela¬ 
tions  are  also  given  by  Caiman  [1978a]  and  Hayashi  [1979],  but 
they  use  a  different  sign  convention  in  defining  their  rotary- 
component  cross-coherence  spectrum,  producing  a  phase  of 
opposite  sign  to  (25). 

For  a  single  wave  of  vertical  wavenumber  m  ,  in  a  given 
sequence  of  profiles,  then  $±(/?!|)  =  2t,  where  r  is  the 
azimuthal  alignment  of  the  wave's  semimajor  ellipse  axis  as 
given  by  the  Stokes  parameters  [e.g.,  Kraus,  1966;  Vincent  and 
Fritts,  1987;  Eckermann  and  Vincent,  1989],  This  was  noted  by 
Gonella  [1972],  who  then  referred  to  |S,(m)|  as  the  stability 
(coherence)  of  the  ellipse  orientation  (see  also  Mooers  [1973]). 
Thus  the  rotary  method  can  be  used  to  investigate  horizontal 
as  well  as  vertical  wave  propagation  directions  and  many  other 
aspects  of  the  wave  field. 

However,  note  once  again  that  the  rotary-component  cross 
coherence  is  not  a  sensitive  measure  of  the  presence  or  ab¬ 
sence  of  a  coherent  wave:  for  example,  a  circularly  polarized 
wave  (P  —  D  =  0)  in  the  presence  of  some  noise  variance 
always  yields  |S*(m)|  =  0,  since  either  the  clockwise  or  the 
anticlockwise  rotary  component  vanishes  in  this  case.  To  iden¬ 
tify  coherent  waves  using  the  rotary-vector  method,  the  spec¬ 
tral  degree  of  polarization  (19)  can  be  calculated  from  the 
rotary-component  spectra  according  to  [Hayashi,  1979] 


1  —  d2(m ) 


4ACW(w)CIT(w) 
[ACW(m)  +  CW(m}f 


[1  -  [S,(m)|2]. 


(26) 


Summary  and  Conclusions 

Continuing  improvements  in  the  accuracy  and  resolution  of 
data  from  atmospheric  profiling  instruments  holds  the  promise 
of  a  wealth  of  new  data  on  atmospheric  gravity  waves  [e.g., 
Hoppe  and  Fritts,  1995;  Gardner  et  al.,  1995;  Hamilton  and 
Vincent,  1995].  Parallel  theoretical  developments  now  mean 
that  spectral  analysis  of  these  data  will  be  a  key  tool  in  the 
appraisal  of  various  theories  of  the  wave  field  [e.g.,  Fritts  and 
Hoppe,  1995;  Gardner  et  al.,  1995].  Some  theoretical  models 
suggest  that  these  wave  spectra  result  from  coherent  large- 
amplitude  waves  with  modulated  intrinsic  parameters  [e.g., 
Sato  and  Yamada ,  1994;  Alexander,  1996;  Warner  and  McIntyre, 
1996],  while  others  conclude  that  they  are  produced  by  an 
incoherent  spectrum  of  many  randomly  interacting  waves  [e.g., 
Hines,  1991;  Medvedev  and  Klaassen,  1995],  In  assessing  these 
and  other  issues  through  data  analysis,  the  coherence,  polar¬ 
ization,  and  energetics  of  the  wave  field  are  pertinent  and  can 
be  quantified  using  a  number  of  standard  analysis  techniques. 
It  is  important  that  the  relationships  among  these  techniques 
are  clear  so  that  firm  conclusions  can  be  drawn. 

By  using  expressions  for  the  power  spectral  densities  of  the 
four  Stokes  parameters  [ Eckermann  and  Vincent,  1989],  rela¬ 
tionships  among  the  autospectra,  cross  spectra,  and  rotary 


spectra  of  wave-induced  horizontal-velocity  perturbations  have 
been  presented.  Since  the  Stokes  parameters  provide  informa¬ 
tion  about  the  polarization  characteristics  of  any  coherent 
gravity  waves  in  the  data,  these  relationships  aid  interpretation 
of  the  data  provided  by  cross-spectral  and  rotary-spectral  anal¬ 
yses.  The  techniques  are  complementary  but  have  different 
strengths  and  weaknesses.  In  particular,  it  was  demonstrated 
here  that  cross-coherence  measures  between  components  of  a 
single-vector  series  are,  in  general,  unreliable  indicators  of  the 
presence  of  a  coherent  wave  in  a  given  data  set:  a  spectral 
“degree  of  polarization”  measure  was  presented  for  future 
studies  using  any  of  the  above  methods. 

There  are  many  more,  potentially  useful,  spectral  quantities 
which  can  provide  further  information  on  the  wave  field  [e.g., 
Mooers,  1973;  Caiman,  1978a,  b;  Hayashi,  1979],  The  informa¬ 
tion  provided  by  these  more  derived  measures  should  also  be 
interpretible  through  the  relationships  among  the  basic  hodo- 
graphic,  cross-spectral,  and  rotary-spectral  quantities  discussed 
both  here  and  elsewhere. 
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